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Abstract 

In this work, we consider the singular integrals of Cauchy type of the forms 

j{f^x)^^^IE^ C — dt, -i<x<i. (1) 

and 



$(/,z) = -^^ / ^=%f- -dt, zi[-l,l]. (2) 

TT J-1 %/l - t^{t - Z) 

which are understood as Cauchy principal value integrals. Quadrature formulas (QFs) for 
singular integrals (Sis) (1) and (2) are of the forms 

N 

J{f,x)=Y,Ak{x)f{tk)+RN{f,x), -l<x<l. (3) 

fe=0 

and 

N 

= + zi[-l,l]. (4) 

fc=0 

where z is complex variable with |i?e(z)| > 1. With the help of linear spline interpolation, we 
have proved the rate of convergence of the errors of QFs (3) and (4) for different classes (i.e. 
iJ"([-l,l],iv:),C""'"[-l,l],T^''[-l,l]) of density function f{t). It is shown that approxi- 
mation by spline possesses more advantages than other kinds of approximation: it requires 
the minimum smoothness of density function f{x) to get good order of decreasing errors. 



1 Introduction 

The importance of singular integrals (Sis) of the form (1) and (2) and their numerical 
solution are given in many researchers work ([l]-[3]) and literatures cited therein. Many 
of them are based on the approximation of density function f{t) with Chebyshev polyno- 
mials. 

Note that in (2), the function ^/ — 1 is understood as a single- valued branch in the 
plane of complex variable with cut along the interval [-1,1] such that a/^^ — 1 = z + 0{z^^) 

^This paper is dedicated to the memory of Professor Israilov M.I. (1934-2010). It was published in collection 
of papers entitled "Differential Equations and Inverse Problems, Press FAN, Tashkent, 1986, pp. 236-258 (in 
Russian) 
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for the large ^ > 0. In the future, under W — arcsinz only a branch of the function for 
which |i?e(VF)| < | will be understood. 

In this paper, we construct efficient quadrature formulas (QFs) for Sis (1) and (2) using 
linear spline interpolation. Obtained QFs provide uniform convergence for any singular 
point X e (—1, 1) and any z ^ [—1, 1]. 



2 Construction of the quadrature formula 



In order to write the exact form of coefficients of the quadrature formula (3) and (4), we 
introduce the following notations 



Gk{x) 



tkVT 



x^ 



irhi 



arcsintk+i — arcsintk), 



■ ztk+i — 1 .ztk 
arcsm arcsm- 



1 



^^(^) = ^f . . 

TT/lfe \ Z — tk+1 Z — tk 

II x ^tk, then the coefficients of QFs (3) are computed by the formulas 

ti — X , 



Ao{x) 
Ak{x) 



-lnGi{x) — — x'^go, 



tk+1 - x^^ Gk+i{x) ^ X - tk-i j^^ Gk{x) 



nhk Gk{x) nhk-i Gk-i(x) 
-Vl-x'^{gk - 9k-i): k = l,...,N-l 



> . 



An{x) = 



tN- 



X 



nh 



N-l 



lnGN-i(x) + Vl - x^qn-i- 



(5) 



As (^^(±1) = 1 for all k, then Ak{±l) = for A; = 0, ...,N. These correspond with the 
fact in [6] that J{f,x)\x=±i = is independent from the value of /(±1) . 

If X coincides with the nodes tk, {k = 1, ...,N — 1), then the coefficients Aj{tk),j ^ 
k — l,k,k + l are computed by (5). If k in (5) is replaced by A; — 1 and k + 1 and x = tk is 
put, then coefficients Ak-i{tk) and Ak+i{tk) are again computed respectively by (5) and 
for Ak{tk) we have 



Ak{tk) 



1 j^^Gk+i(tk) 

TT Gk-l{tk) 



- \n-tl{gk-9k-i), k^l,...,N-l. 



(6) 
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Coefficients of tfie QFs (4) fiave the form 



Bo{z) = {z- h)Fo{z) + Vz^ - Igo: 

Bk{z) ^ {z- tk+i)Fk{z) -{z- tk-i)Fk-i{z) 

+Vz^-l{gk-gk-i), k^l,...,N-l 



> . 



Bn{z) = {tN-i - z)Fn-i{z) - Vz^ - lgN-1- 



(7) 



Let us derive coefficients of QFs wliicli are given by (5) and (7). As we know tlie linear 
spline SN{t) interpolating the given function / on the grid A: — l = to<^i<---< 
^AT-i < = 1 for t e [tj^tj+i] has the form 



Ilk L 



(8) 



Replacing f{t) in (1) with SN{t) we have 



J{Sn,x) = 



N-l 

El 



k + 1 



n — ' hk 



{tk+i - t)f{tk) + {t- tk)f{tk+l) 



tk 



VT^{t - x) 



dt. 



x^ 



TT 



ho J VT^{t 



/(«+ ' 



X] 



{t — tisf-i)dt 
hN-i J Vl - t^{t - x) 

tN-l 



N-2 



k=l 



{t - tk)dt 



*fc+2 



{tk+2 - t)dt 

y/T^(t-x) ' J VT^(t-x) 
*fc+i 



+ 



fitk+i). 



(9) 



Introducing notations 



*fc+i *fc+i 



tdt 



tk 



VT^it-xY 



and using easy checking formulas 

xdt 



I 



tdt 



, 1 + VT^^y/T^ + Xt ^ 
In , , h G, 

= -ln{VT^ + VT^) + c, 
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obviously we have 

J{k, x) 



1 



-In 



{t - + Vl-x^VT^ + xt) 



*— *fc+i 
t=tk 



writing 

{t - + Vl - a;Vl - + xt) = (Vl - x^ + Vl - t'^){tVl - x^ - xVl-t^), 
we obtain 



J{k, x) — 



1-x' 



t=tfc+i 



t=tk 



1 , Gk+i 
rln- 



1-x^ GkixY 



Next, having the foUowing relation 

Ji{k, x) = xJ{k, x) + nhkgk- 

Eq. (9) can be rewritten as 



(10) 



(11) 



J{Sn,x) 



VI -x^ r 1 



TT 



hn - 



tiJ{0,x)-Ji{0,x) /(-I) 



N-l 



— \ji{N -l,x)- h-iJi{N - 1, x)] /(I) 



+ Y{^k+iJ{k,x) - Ji{k - l,x)^ 

k=i L 

+ 7^ (Jiik -l,x)- tk-iJ{k - 1, x)) 



f{tk) 



(12) 



Substituting (10) and (11) into (12) and simplifying the expressions, we arrive at (5) for 
finding the coefficients Ak{x) of the QFs (3). Furthermore, we can derive the Eq. (6) 
from (12) and (lO)-(ll) as follows 



Akitk) 



VI -x^ f 1 
h 



-(tk+iJ(k,tk+i) — Ji(k,tk)J 
(ji(A; - 1, tk) - tk-iJ{k - l,tk)) I 



h 



TT 



[J {k, tk) - J{k-1, tk) - n{gk - Qk-i)] 



1 ^^G^ - \/l - tl{9k - 5.-i), k = 1, ...,N - 1. 

TT Lrk-l[tk) 



(13) 
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In order to derive the coefficients Bk of tfie QFs (4) we use tlie combination of tfie integrals 

tk+i *fc+i 

r{k,z)= I , Jl{Kz)= [ -^=J£: -. (14) 

^ ' J VT^{t-z) ^ J VT^{t-z) ^ ' 

First J*{k,z) is computed at z = x, where x is any number such that \x\ > 1. Then 
continuing analytical function J*{k,x) along the intervals (— oo, —1), (1, oo) on the plane 
of complex variable z with cut along the interval [—1,1], we obtain 



T*fU \ ^ . zt-1 

J [k, z) — —^=^= arcsm 



" Fk(z). (15) 

- / I 

t=tk 



z-t 
For J* 

Jl{k, z) = zJ*{k, z) + Qkuhk. (16) 

Replacing x into z in (12) and using (14)-(16), we obtain the coefficients Bk of the QFs 
(4) 

3 Estimation of errors 

Let us introduce the following classes of functions: 

1. if"([— 1, 1], i^) is a class function satisfying Holder condition on the interval with 
the index a and constant K. 

2. c™'"[-i,i] = [f{t) : /(-) e m{[-iMKm)] 

3. CC^'"[-1,1] = [f{t) : f{t) e C[-l,l] and f{t) G }. 

4. iy[-l, 1] = < fit) : f^''~^\t) is absolutely continuous and ess sup |/(')| = Mr 

y iti<i 

5. CW'^[-l,l] = [fit) : fit) e C[-l,l] and fit) e ^^[t^, tfe+i] }■ 

Everywhere we use the notation ||/||c = ||/(^)||c[-i,i] as a norm of the function. Note 
that Mr = for any f{t) e C^[-l, 1]. 

Now we prove the following theorems with respect to QFs (3) and (4). 

Theorem 1 Let fit) he a function belonging to one of the classes of functions 1, 1], 
CC]^ or CWlX-l, 1]. Then for the errors of QFs (3) the estimations 

, , „ , , , , , , LuN 
\\RN{f.x)\\c<L^, 

are true for all x G (—1,1), where L and (3 are given in the Table 1. 
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Table 1: For QFs (3) 



Classes of functions 


8 


L 


W'[-l,l] 


1 


47M1 / ^V2\ 




l + a 


TT 27/niV j 




2 


TT \ ' jlnN j 



Remark 1: In the case of uniform grids, 7 = 2. 

Theorem 2 Let f{t) satisfy the conditions of Theorem 1. Then the errors of QFs (4) 
are 

max|iJ^(/,.)|<L-^. L' = ^L'^[^)' 

where L, Li and j3 are given in the Table 2. 

Table 2: For QFs (4) 



Classes of functions 


p 


L 


Li 


W\-l,l] 


1 


47Afi / ^^2 \ 
TT 2jlnN j 


Mi7 
27r 


CCi"[-l,l] 


l + a 


TT \^ 27/nAr j 


K17I+" 
Air 


CWl[-l,l] 


2 


72M2 / ^ 7rV2 \ 
TT \ jlnN j 


M272 
Stt 



Next theorem is again related to QFs (3) but in different classes of functions: 
Theorem 3 Let f{t) be a function belonging to one of the classes of functions 

H-{[-l,l],K), W'[-l,l], CC]^" orCWli-1,1]. 
Then the error terms of QFs (3) satisfy the following estimations 

InN 

\\RN{f,x)\\c < 

for all X E (—1,1), where L2 and (3 are given in the Table 3. 



Remark 2: Note that the main terms of L2 in the Theorem 3 for the last three classes of 
functions is twice less than the main terms of L in the Theorem 1. 
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Table 3: For QFs (3) 



Classes of functions 



/3 



CCi"[-l,l] 
CWll-1,1] 



a 



1 + a 



IT 



27M1 



vr 

TT 

27r I 



1 1 + 



1 + 



127r \ 
-flnN J 

127r 

ylnN 



1 + 



7rv247r 



7/niV 



In estimation of the error of QFs (3) we use the idea of [4] (see also [3]) and the 
following Lemmas. 

Lemma 1 Let SN{t) be linear spline (8) interpolating f{t) on the grid A, and let t e 
Then for the estimate of error rN{f,t) — <S'jv(i) — f{t), we obtain 

\\rN{f,t)\\c = r%{hk), 
where r*j^{hk) is defined by the Table (4) 

Table 4: Error of the linear spline (8) 



Classes of functions 










\hk 






cwl[-i,i] 





Lemma 1 is proved as Theorem 2.1 which is shown in [4]. 

Lemma 2 Let SNit) be linear spline defined by (8). Then 

rN{f.t)eH\[-lMk), 

where K is given in Table (5). 
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Table 5: Error of the linear spline (8) 



Classes of functions 


K 




2Mi 






CWli-1,1] 


M2hk 



Proof of the Lemma 2. Consider three cases: 

{a)t,t' e [tk,tk+i], {b)\t-t'\>hk {c)t e [tk-i,tk], t' e [tk,tk+i], \t-t'\<hk 
I. Let f{t) e 1, 1]. Then in the case (a), using representation (8) we have 



rN{f;t)-rN{f;t') 



1 

hk 



1 

hk 

ik+l 



{t-W{tk+i)-f{h)]-hk[f{t)-f{t')] 



{t - 1') J f{s)ds -hkj ns)ds 



< 2Mi\t-t'\ 



In the case (6), in accordance with Lemma 1, we get 

riv(/;t) -r^(/;t')| < \rN{f]t)\ + |rAf(/;t')| < hkMi < M^\t-t'\. 
Using the case (a), in the case (c) we obtain 

rjv(/; t) - rjv(/; t') I < |rjv(/; t) - TNif; tk) | + |r^(/; 4) - r^if; t') 
<2Mi\tk-t\+2Mi\t-t'\ ^2Mi\t-t'\. 

II. Now let f{t) e CCi"[-l, 1]. In the case (a) 

1 



riv(/;t)-r^(/;0 



< 



hk 
1 

hk 



{t-t')[f{tk+i)-f{tk))-hk[f{t)-m] 

{t - t'){tk+i - tk)f'{9i) - hk{t - t')f{92) 



= {t-t') /(^i)-/(^2) <Ki\t-t'\\ei-e2r 

< K^hl\t-t'\. 



In the case (6), due to Lemma 1, we have 

r7v(/;t)-r^(/;0 <K^hl\t-t'\ 
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It is obvious in the case (c) that 

rN{f;t)-Mf;t')\<K,ht\t-t'\. 
III. Let f{t) e CWl[-l, 1]. In case (a), we have 



rN{f;t)-rj,{f;t') 



t-t' 
t-t' 



no,) - m) 

I ns)ds 
01 



< M-2hk\t-t'\. 



The cases (6) and (c) are proved in a similar way as the case (a). So that the proof of the 
Lemma 2 follows from the above obtained errors. 

Now it is easy to prove the following lemma. 



Lemma 3 Let SN{t) be linear spline defined by (8) and f{t) e 1, 1],K). Then 



Prove of the Theorem 1. Since 



1 

/ 



dt 



VT^^{t - x) 



= 0, 



the reminder term of QFs (3) can be represented as 



TT 



-1 



Vl-x^{t-x) 



(17) 



For definiteness, let us prove the Theorem 1 in case < x < 1 (the case — 1 < x < is 
considered analogically). Fixing the number < 5^ < ^ and dividing the integral in (17) 
into three parts to yield 



RN{f,x) = 



X— 5jv 1 



TT 



-1 x—6n x+Sm 



rN{f,t) - TnU^x) 
VI -a;2(t - x) 



dt 



x^ 



TT 



[Jl + J2 + h) ■ 



(18) 
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First assume that Sn < ^ — x. Then due to (10), for Ji we have 

x—Sn 

Jl 



A\rN{f,x)\\c 
2 



dt 



It is not hard to show that 



rjv(/, x)\\cln 



VT^^{t - x) 



t=x—S]\ 



t=x—S\ 



1 



(l-x''- v^r^^Vl -{x- Sn^^ . 



^ X + — [1 — X 
On V 

This is a function of x which strictly decreases on [O, ^] and strictly increases on l] , 
and 



i + ^r^ 



N 



Hence, 



l^il < 



x^ 



\rN{f, x)\\cln—. 

On 



(19) 



For J2, we use Lemma 2 



IJ2I < K 



dt 



x—Sf, 



i^[arcsin(x + 5^) — arcsin(x — 5^)]- 



Let 



(fi2{x, S]sf) — arcsin(x + Sn) — arcsin{x — 5n)- 



Since < x < 1 — 5n and by assumption 5n < i — x, derivative of (fi2{x, 5n) is positive 
and 

ip2{x, Sn) < </22(l - Sn, Sn) = arcsin2A/5Ar(l - Sn)- 
From this and the known inequality arcsino; <|a,(0<Q;<|)it follows that 

^2ix, Sn) < n>/S^. 

Hence 

\J2\<K^T^/S^ (20) 
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For J3, we have 

|J3| = 2||r^(/,x)||c 



dt 



Mf,x)\\c{-l)ln 



t=X+dM 



We may show that the function 



<P3{x, 5n) 



= -(pi{x, -Sn), 



t=X+5N 



strictly increases on [0, 1 — ^at] and strictly decreases from — ^.^ to 1. So that 



IJ3I < 



X' 



\rN{f:X)\\cln—. 

Or 



On 



It follows from the errors of (19)- (21) and (24) that 



\RNif, x)\\c<-\ Mf, x)\\cln^ + Ky/S^. 



Now consider the case Sn > ^ — x. Write 



-Rjv(/, x) 



x^ 



TT 



' x— <5jv 1 
+ 

-1 x—Sn 



rN{f,t) -rN{f,x) 
\J\ — x'^it — x) 



dt 



X' 



TT 



( jr + J*) 



Integral Jj* is estimated as Ji. Due to Lemma 2 

1 



\j;\<K 



d.t 



x—S]\f 



~ K arcsin a/I — (x — S^y. 



Since 0<x<l, 5jq > 1 — x and due to the inequality 

l-{x- d^f = {1-x + 6n){1 + x-8n)< 45jv- 

we obtain 



IJ2I < -^arcsin 2 a/^jv < Kt:\/5 



'N- 



(21) 



(22) 



(23) 



Substituting the errors of J* and into (23), we arrive at estimation (22). 

In order to determine the errors of estimation for every classes of functions in Theorem 
1, we use the results of Lemma 1 and 2 and set = Viz: 
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I. Let f{t) e 1]. Then 

\\RN{f,x)\\c <-M,hln^ + 2M,^/S^ 

TT On 

_ 4Mi7 / 7r\/2 \ In 

II. lif{t) e CCi"[- 1,1], then 

I |i?jv(/, x)\\c< — In ^ + 2K^h''^ 

_ 2i^i7i+" / 7r\/2 \hiN 

~ ^ y ^ 27 In AT J A^i+"' 

III. If f{t) e CVFl[-l, 1], then 

IIP (f ^11 ^ ^2 7% 2 2M27 



TT I 7lnAr / ATS 



Theorem 1 is proved. 



Proof of the Theorem 2 is carried out by the famous scheme of the formula 
Sokhotskii-Plemergh (see [5]), principle maximum module for analytical function and 
results of Theorem 1 and Lemma 1. 

Proof of the Theorem 3. Let the remainder term of QFs (3) be divided into three 
parts 

(x— 5jv X+&JSI 1 \ 
J ^ J h Vi-x-(i-x) 

= ^^^IZZ + + J3) . (24) 

In the proof of Theorem 1, we have already seen that the case ^at < 1 — is adequate 
for the estimations of Ji and J3 i.e. 

4 



IJ1I + IJ3I < -Vl-x^\\rN{f,x)\\cln—. (25) 
TT Oat 



For the estimation of J2, we consider the function 

x+ve 



\t — x\ 



T{x,e,a) = Vl^ I ^^-j^=-dt, (26) 
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where 0<a-<l, 1- x > s, 0<£<^. 



It is obvious that 



T{x,e,a) < 28" 



1/2 



-.dy + 



VI - x^y"-^ 



V"! - (x + £y)2 J ^\ - {x + eyy 

1/2 

= 2£'^[Ti(x,£,a)+r2(a:,£,a)]. 
Let fco = [1/^]) since 1 — x > s, then for some k,l < k < ko the inequahty 

A;£ < 1 - X < (A; + l)e, 

takes place. From this and < y < - it follows that 



1 - (x + £y)2 (1 - X - £y)(l + X + £y) ~ l-x-s/2 ~ A; - 1/2 
for all k > 1. Hence 

1/2 

1— a 



T^{x,e,a) = 2 j y-^dy^'^. 



Furthermore 



n{x,e,a)^^l -= 



VT^ 



1/2 



\^T^^JxTey) 



dy. 



Let £ < 1 - X < 2£. Then 



1-x^ 1-x 2 
< < 



1 — {x + eyy 1 — X — I 1 — 



and therefore 



T,{x, e, a) = 2'-^V2 J {1 - yf'^dy = 2^"^ 
If for some k <2 the inequahty 



1/2 



A;£ < 1 - x < (A; + 1)£ 
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takes place, then 

l-x^ 1-x k + 1 
< < < 3 

1 — {x + syY ~1 — X — ey~k — 1" ' 

therefore 



T2 {x, e,a)< 2^~''VS ■ \ < 2^'''. (30) 



Prom (27)-(30) it follows that 



r(x,£,a) < 2^-^ ( 2 + - ) e'^. (31) 



a 



Now for J2 we have 



j2<-vT^ / ' ^-r(x,£,a), (32) 



a;— (5 



where cr = a for the class if"([— 1, 1], X) and cr = 1 for rest classes of functions. Assuming 
o — and from the Lemmas 2, 3 and inequalities (25) and (31)-(32) we get the assertion 
of the Theorem 3. 
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